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We present a new search for H°V production, where H° is a scalar Higgs boson decaying into 
bb with branching ratio (3, and V is a Z° boson decaying into e + e _ , p + p~ , or vv. This search is 
then combined with previous searches for H°V where V is a boson or a hadronically decaying 
Z°. The data sample consists of 106 ± 4 pb -1 of pp collisions at y/s =1.8 TeV accumulated by the 
Collider Detector at Fermilab. Observing no evidence of a signal, we set 95% Bayesian credibility 
level upper limits on a(pp —> H°V) x (3. For H° masses of 90, 110 and 130 GeV/c 2 , the limits are 
7.8, 7.2, and 6.6 pb respectively. 

PACS numbers: 14.70.-e, 13.85.Qk, 13.85.Ni 


A key component of the standard model (SM) is spon¬ 
taneous electroweak symmetry breaking, which gives rise 
to the mass of all fermions and the W ± and Z° gauge 
bosons. This process leads to the existence of a neu¬ 
tral scalar particle, the Higgs boson (if 0 ), whose mass is 
unspecified in the SM, but whose couplings to all other 
particles of known mass are fully specified at tree level. 
The Higgs boson has not been directly observed, but its 
expected contribution to loop corrections for many SM 
observables has allowed an inferred mass of Mh = 126 l 4 g 
GeV/c 2 from precision electroweak measurements 0. In 
addition, direct searches at LEP 2 have excluded, with a 
95% confidence level, a SM Higgs boson with Mh < 114.4 
GeV/c 2 J5]. The relatively low H° mass favored within 
the SM framework implies the possibility of its direct ob¬ 
servation at the Tevatron in Run II, where searches have 
begun by both the D0 Q and CDF collaborations. Here 
we report on direct searches using data accumulated by 
the Collider Detector at Fermilab (CDF) between Febru¬ 
ary 1992 and July 1995 (Run I) for a total integrated 
luminosity of 106 ± 4 pb -1 . 

At the Tevatron, the SM Higgs boson is produced from 
both gluon-gluon and quark-antiquark initial states 0 . 
Although the dominant production mechanism is gg —> 
H°, production in association with vector bosons (qq' —> 
qq — > H°Z°) provides the most sensitive chan¬ 
nels for Higgs boson searches at the Tevatron if Mh < 
140 GeV/c 2 , because one can obtain significant back¬ 
ground rejection from the additional highly energetic ob¬ 
jects in the event coming from the vector boson decays. 
The predicted cross section, a VH o, for VH° production 
from pp collisions at yfs = 1.8 TeV varies between 0.50 
and 0.15pb for H° masses between 90 and 130GeV/c 2 , 
with the ratio <j\y±H° /&z°h 0 ~ 1 . 6 . 

We have previously reported the results of searches in 
the WH° — > £vbb (£ = e or p) and VH° — > qq 1 bb chan¬ 
nels Si Here we add the searches for Z°H° production 
using the decay channels £ + £~ bb and vv bb. Finding no 
evidence for Higgs boson production using these decay 
modes, we set limits on the production cross section as a 
function of mass, and combine our results with the previ¬ 
ous VH° cross section limits. These limits represent the 
final CDF cross section limits for Higgs boson production 


in association with a vector boson from the Run I data. 

The CDF detector is described in Ref.jV], and the co¬ 
ordinate system and various quantities used throughout 
this paper are defined in Ref. jlj. The momenta of the 
charged leptons are measured with the central tracking 
chamber in a 1.4 T superconducting solenoidal magnet. 
Electromagnetic and hadronic calorimeters surrounding 
the tracking chambers are used to identify electrons and 
jets and measure their energies. Muons are identified 
with drift chambers located outside the calorimeters. 
The silicon vertex detector (SVX) is the innermost de¬ 
tector used for precise tracking in the plane transverse to 
the beam @]. 

In the analyses reported here, two algorithms using 
tracks measured with the SVX are applied to identify jets 
originating from heavy flavor quarks (b and c). The first 
reconstructs a secondary vertex (a vertex displaced from 
the primary interaction vertex) produced by the heavy 
flavor decays and measures the transverse decay length 
(SVX tag). The resolution of the transverse decay length 
of the secondary vertex is typically of order 150 pm. The 
second algorithm uses the impact parameter of the tracks 
in the jet (the closest distance of the track to the primary 
vertex in the transverse plane) to calculate a probability 
that the jet is not from heavy flavor (JPB tag). The 
details of these tagging algorithms are given in Ref. E3- 

For the details of the analyses previously published we 
refer to those publications [g, |6|, and list the results here 
when appropriate, as they are used in the combined cross 
section limits. We now describe the two new channels, 
l + £~ bb and vv bb. Events for the £ + £~ bb channel analy¬ 
sis are required to pass a high-P t lepton trigger and must 
contain two high-Pr [s], oppositely charged leptons (e or 
p) that are isolated from nearby tracks and calorimeter 
activity. At least one lepton is required to have Pp > 20 
GeV/c and be in the central detector (\rj\ < 1.0). For the 
second lepton the Pp requirement is relaxed to 10 GeV/c 
and the pseudorapidity range is extended into the plug 
calorimeter, up to \rj\ ~ 2.4. The dilepton invariant mass 
must be in the range 76 < Mu < 106 GeV/c 2 to be 
consistent with the decay of a Z° boson. This require¬ 
ment essentially removes any sensitivity of this analysis 
to Z° — > t + t~ . The event is additionally required to 
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contain two or three high-F t jets ( Et > 15 GeV), at 
least one of which is SVX tagged. A cut on the miss¬ 
ing transverse energy (Jp T < 50 GeV) [§j is also applied, 
with the effect of reducing the tt background by approx¬ 
imately a factor of two, while preserving about 95% of 
the signal. 

The vv bb channel is characterized by two heavy flavor 
jets and large Ifir from the neutrinos. The data sample 
for this channel is derived from an event trigger requiring 
IpT > 35 GeV in addition to event quality cuts |llj. To 
reject W ± and Z° decays to leptons, events containing 
an isolated track with Pt > lOGeV/c are removed from 
the sample. To ensure less susceptibility to the uncer¬ 
tainty in the trigger efficiency at threshold, the analysis 
requires Ifix > 40 GeV. The trigger efficiency is approxi¬ 
mately 60% at this value of Ifix- Additionally the event 
must contain two or three jets with Et > 15 GeV (about 
20% of the ZH° signal contains a third jet). To reject 
QCD multi-jet events where the results from a mis- 
measured jet, the azimuthal angle between the fr and 
the direction of any jet with Et > 8 GeV is required to be 
at least 1.0 radians. In addition, the jets from inclusive 
di-jet production tend to be back-to-back, while jets from 
H° —> bb in ZH° events tend to have a smaller opening 
angle, leading to the requirement that the azimuthal an¬ 
gle between the leading two jets be less than 2.6 radians. 
Approximately 10% of the efficiency from the vv bb selec¬ 
tion is contributed by W^H 0 events where the lepton is 
undetected. 

Events in the uDbb sample are classified as “single- 
tagged” (exactly one SVX tagged jet) or “double-tagged” 
(one SVX tagged jet and a second jet tagged by either 
the SVX or JPB tagging algorithms). The backgrounds 
and efficiencies are calculated separately for these orthog¬ 
onal sets, which are then treated as separate but corre¬ 
lated channels when combined with the other channels. 
This is analogous to what was done in the WH° —> ivbb 
search jg). 

TABLE I: Total selection efficiencies for VH° events in each 
analysis channel used in the combined result, as a function of 
the H° mass, Mh (GeV/c 2 ). Numbers are percentages and 
include the branching ratios of the vector boson (W ± or Z°) 
in a given channel. ST refers to single-tagged events and 
DT to double-tagged events. Uncertainties include systematic 
effects. 


Channel 

VH° 
Mh = 90 

event efficiencies (%) 

Mh = 110 Mh = 130 

i+r bb 

0.14 ±0.03 

0.20 ±0.04 

0.19 ±0.04 

vvbb (ST) 

0.51 ±0.10 

0.63 ±0.13 

0.76 ±0.15 

vvbb (DT) 

0.37 ±0.08 

0.43 ±0.09 

0.51 ±0.10 

ivbb (ST) 

0.59 ±0.15 

0.72 ±0.18 

0.80 ±0.20 

ivbb (DT) 

0.22 ± 0.06 

0.29 ±0.07 

0.30 ± 0.08 

qq'bb 

1.3 ±0.7 

2.2 ± 1.1 

3.1 ± 1.6 


The efficiencies for identifying VH° events with our 


selection criteria are summarized in Table [I] and are de¬ 
termined from a PYTHIA 0 Monte Carlo simulation of 
Higgs boson production via V* —> VH° —> Vbb followed 
by a detector simulation. The Higgs boson is forced to 
decay to bb with a 100% branching ratio. The identifi¬ 
cation efficiencies for single leptons are measured from 
Z° —> i + i~ events in the data and are found to be 
91% for muons and 83% for electrons [l3|. The SVX 
and JPB 6-tagging efficiencies are determined using data 
and Monte Carlo samples with high 6-purity 0- In the 
£ + £~ bb channel, the efficiency for obtaining > 1 SVX tag 
in a signal event is (45 ± 7)%. The double 6-tagging effi¬ 
ciency in the vv bb channel (SVX + SVX or SVX + JPB) 
is (19±4)%, and the single 6-tagging efficiency (one SVX 
tag) is (25 ± 3)%. The total event efficiencies are the 
product of the trigger efficiencies, the kinematic and ge¬ 
ometric acceptances from the selection cuts, the lepton 
identification efficiencies when appropriate, the 6-tagging 
efficiencies, and the V branching ratio relevant for a given 
search channel. The systematic uncertainties in the total 
efficiencies for the £ + l~ bb and vv bb channels are approx¬ 
imately 20%, comprised mostly of the uncertainties in the 
6-tagging efficiency (15%), the modeling of initial and fi¬ 
nal state radiation (7%), lepton identification efficiency 
(7% for the £ + i~ bb channel), trigger efficiency (5% for 
the vv bb channel), integrated luminosity (4%) and the 
energy scale of jets (3%). 

In the £ + £~bb channel the dominant background is Z° 
production in association with a heavy flavor pair ( Z°bb , 
Z°cc), which accounts for approximately 60% of the to¬ 
tal. About 20% comes from Z°+jets events where a jet 
is mistagged due to track mismeasurements, and there 
are smaller contributions from Z°c , Z°b, diboson, and tt. 
All backgrounds are determined using Monte Carlo sim¬ 
ulations except that from Z°+jets which uses the data. 

The vvbb channel background is dominated by QCD jet 
production of bb where the $t results from mismeasured 
jets. To calculate this contribution we first parameterize 
the tagging rate in the t < 40 GeV region of the data as 
a function of jet Et and track multiplicity. By applying 
this parametrization to the jets in the signal region we 
estimate the QCD background to be about 70% of the 
total background in the single-tagged sample and about 
50% in the double-tagged sample. Smaller backgrounds 
include V+ heavy flavor, diboson, and tt, all of which are 
derived from Monte Carlo simulations. 

For each decay channel, Table UH summarizes the to¬ 
tal expected backgrounds, the expectations from stan¬ 
dard model VH° production for Mh = 110 GeV/c 2 
and H° —> 66, and the number of data events observed. 
The dominant background in the qq'bb channel is QCD 
production of bb with additional jets, hereafter abbrevi¬ 
ated as “QCD”. Its normalization is difficult to predict 
and therefore left unconstrained in the analysis. Fur¬ 
ther details of the background calculations are given else¬ 
where 11110 








5 


TABLE II: Predicted numbers of events in each channel from 
all backgrounds (see text), expected number of signal events 
for Mh = 110 GeV/c 2 , and number of events observed. Un¬ 
certainties include systematic effects. There is no reliable pre¬ 
diction for the background in the qq bb channel. 


Channel 

Background 

Signal 

Data 

£ + £~ bb 

3.2 ±0.7 

0.06 ±0.01 

5 

vvbb (ST) 

39 ±4 

0.20 ± 0.04 

40 

vubb (DT) 

3.9 ±0.6 

0.14 ±0.03 

4 

£ubb (ST) 

30 ±5 

0.23 ± 0.06 

36 

£ubb (DT) 

3.0 ±0.6 

0.09 ±0.02 

6 

qq bb 


0.73 ±0.29 

589 


A binned likelihood is used to compare the dijet mass 
spectrum (of the two tagged jets, or the one tagged jet 
and highest-Er untagged jet) in the data to a combina¬ 
tion of expected distributions from the background pro¬ 
cesses and the VH° signal, as a function of H° mass. 
The observed dijet mass spectra for the vv bb and £ + £~ bb 
channels are shown together with the expected back¬ 
ground shapes in Figs. 0 and 0 respectively. 
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FIG. 1: Dijet invariant mass in uvbb candidate events, for 
events with exactly one 6-tagged jet and separately for events 
with two 6-tagged jets. The single 6-tag data includes one 
overflow event. The background shapes shown differ only in 
the inclusion of the predominant background of QCD bb pro¬ 
duction. The signal shape shown (dashed line) is for a SM 
Higgs mass of 110 GeV/c 2 and a normalization of 50 times 
the expected rate. 


Since no signal is observed, we calculate upper limits 
on VH° production using a Bayesian procedure. For each 
channel, a posterior density is obtained by multiplying 
the likelihood function for that channel with prior den¬ 



FIG. 2: Dijet invariant mass in £ + £~ bb candidate events. 
At least one jet is required to be 6-tagged by the SVX al¬ 
gorithm. The background shapes shown differ only in the 
inclusion of the predominant background of Z°±heavy flavor. 
The signal shape shown (dashed line) is for a SM Higgs mass 
of 110 GeV/c 2 and a normalization of 50 times the expected 
rate. 


sities for all the parameters in the likelihood: integrated 
luminosity, background normalizations, signal efficiency, 
and the product a ' VH0 = (Jvh° x (3 of the signal cross sec¬ 
tion (Tyro by the branching ratio (3 for H° —> bb. With 
two exceptions, these priors are truncated Gaussian den¬ 
sities constraining a given parameter to its expected value 
within its uncertainty. The exceptions are <j' V ro and 
the QCD background normalization in the qq' bb chan¬ 
nel. Since nothing is presumed known a priori about 
these parameters, they are assigned uniform priors. The 
posterior density is then integrated over all parameters 
except cr' VH0 , and a 95% credibility level (C.L.) upper 
limit on <?' VH0 is obtained by calculating the 95 th per¬ 
centile of the resulting distribution. When combining 
channels, the same procedure is applied to the product 
of their likelihoods. Correlations in the total efficiencies 
are taken into account by identifying common parame¬ 
ters such as the 6-tagging efficiency and some kinematical 
efficiencies. Each of these common parameters is then as¬ 
signed a single prior. 

Upper limits on ay#o x (3 in each channel and in all 
channels combined are summarized in Table HTTl as a func¬ 
tion of H° mass. 

These results are also plotted in Figure 0 The stan¬ 
dard model prediction is about 30 times smaller than the 
measured 95% C.L. upper limits. For the Ivbb and qq' bb 
channels, the limits reported here are slightly different 
from those previously published @, @]; this is mainly due 
to our improved understanding of the 6-tagging efficiency 
(ltij. Table HTTl also shows expected upper limits under the 
assumption of zero signal. These expectations are calcu¬ 
lated over an ensemble of experiments similar to this one, 
but where the background normalizations are fluctuated 
around their expected values by their uncertainties. We 
note that the observed combined limits are driven by the 
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TABLE III: The 95% credibility level upper limits on a(pp —> 
VH°) x P where /? = BR (H° -> 66), for each of the search 
channels and their combination, as a function of H° mass, 
Mh (GeV/c 2 ). Also shown are the expected limits under the 
assumption of no H° signal. ST designates the single 6-tagged 
subsample and DT the double 6-tagged subsample. 


Channel 

Measured (expected) upper limits (pb) 
M h = 90 M h = 110 M H = 130 

£+e~ bb 

55.6 (36) 

31.8 (24) 

23.8 (25) 

vvbb (ST) 

20.8 (30) 

20.8 (21) 

18.4 (17) 

vvbb (DT) 

10.4 (17) 

9.2 (14) 

8.0 (12) 

vvbb (ST+DT) 

7.6 (13) 

7.8 (11) 

7.4 (8.8) 

ivbb (ST) 

30.0 (18) 

29.4 (15) 

27.6 (12) 

ivbb (DT) 

31.0 (24) 

26.6 (19) 

24.2 (18) 

ivbb (ST+DT) 

23.2 (13) 

22.6 (11) 

21.6 (9.0) 

qq' bb 

38.2 (77) 

21.2 (43) 

17.8 (29) 

All combined 

7.8 (7.1) 

7.2 (5.7) 

6.6 (4.7) 



FIG. 3: Summary of all Run I CDF 95% credibility level 
upper limits on a(pp —> VH°) ■ /3. The lines for the vvbb 
and ivbb channels represent the combined limits from the 
single 6-tagged and double 6-tagged subsamples. Shown for 
comparison is the standard model prediction, and the region 
excluded by the LEP experiments. 


vv bb channel, as a result of an observed downward fluc¬ 
tuation in the dijet invariant mass region of interest in 
this channel. However, the expectation is for the vv bb 
and Iv bb channels to have comparable sensitivity. 

In conclusion, we have searched for Z°H° production 
using the i + i~ and vv decay channels of the Z° and 


produced limits on VH° production using these chan¬ 
nels. We combined these limits with those previously 
published using other decay channels of the vector bosons 
to obtain final CDF Run I 95% C.L. limits on ctvh° x P 
ranging from 7.8 pb to 6.6 pb for H° masses of 90 GeV/c 2 
to 130 GeV/c 2 . These limits additionally apply to any 
scalar particle decaying to bb that is produced in associ¬ 
ation with a vector boson. These results and the com¬ 
bination methodology establish the foundation for our 
searches in the Tevatron Run II data at y/s = 1.96 TeV 
which are exploiting more search channels, an improved 
detector, and more advanced analysis techniques Il4l . 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their contributions. 
This work was supported by the U.S. Department of En¬ 
ergy and National Science Foundation; the Italian Isti- 
tuto Nazionale di Fisica Nucleare; the Ministry of Sci¬ 
ence, Culture, and Education of Japan; the Natural Sci¬ 
ences and Engineering Research Council of Canada; the 
National Science Council of the Republic of China; and 
the A. P. Sloan Foundation. 


[1] The LEP Electroweak Working Group, 
http://lepewwg.web.cern.ch/LEPEWWG/ V.M. Abazov 
et al, Nature 429, 638 (2004). 

[2] R. Barate et al., Phys. Lett. B 565, 61 (2003). 

[3] V.M. Abazov et al., hep-ex/0410062, submitted to Phys. 
Rev. Lett.. 

[4] A. Stange, W. Marciano and S. Willcnbrock, Phys. Rev. 
D49, 1354 (1994). 

[5] F. Abe et al, Phys. Rev. Lett. 79, 3819 (1997). 

[6] F. Abe et al, Phys. Rev. Lett. 81, 5748 (1998). 

[7] F. Abe et al, Nucl. Instr. Meth. Phys. Res. A 271, 387 
(1988). 

[8] In the CDF coordinate system, 9 and <j> are the polar 
and azimuthal angles, respectively, with respect to the 
proton beam direction (z axis). The pseudorapidity p is 
defined as — lntan(0/2). The transverse momentum of a 
particle is Pt = Psin9. The analogous quantity using 
calorimeter energies, defined as Et = E sin 9, is called 
transverse energy. The missing transverse energy, Ifir, is 
defined as — ^ E\. hi, where hi is the unit vector in the 
transverse plane pointing from the interaction point to 
the energy deposition in calorimeter cell i. 

[91 D. Amidei et al., Nucl. Instr. Meth. Phys. Res. A 350, 
73 (1994). 

[10] T. Affolder et al, Phys. Rev. D64, 032002 (2001). 

[11] T. Affolder et al, Phys. Rev. Lett. 88, 041801 (2002). 

[12] T. Siostrand, Comput. Phys. Commun. 82, 74 (1994). 
We use PYTHIA V5.7. 

[13] F. Abe et al, Phys. Rev. Lett. 80, 2779 (1998). 

[14] L. Babukhadia et al. (CDF and D0 Higgs Sensitiv¬ 
ity Working Group), FERMILAB-PUB-03-320-E (2003); 
C. Neu (CDF collaboration), FERMILAB-CONF-04- 
027-E (2004). 



















